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The current status of big bang nucleosynthesis is reviewed with an emphasis on the 
comparison between the observational determination of the light element abundances of 
D, 3 He, 4 He and 7 Li and the predictions from theory. In particular, we present new 
analyses for 4 He and 7 Li. Implications for physics beyond the standard model are also 
discussed. Limits on the effective number of neutrino flavors are also updated. 

1. Introduction 

The overall status of big bang nucleosynthesis is determined by the comparison of 
the rather slowly changing theoretical predictions of the light element abundances 
and the sometimes quickly changing observationally determined abundances. The 
observed elements, D, 3 He, 4 He, 7 Li, have abundances relative to hydrogen which 
span nearly nine orders of magnitude. By and large, these observations are consis- 
tent with the theoretical predictions and play a key role in determining the con- 
sistency of what we refer to as the standard big bang model and its extrapolation 
to time scales on the order of one second. Here, we will review the status of this 
consistency. We begin by briefly outlining the key sequence of events in the early 
Universe which leads to the formation of the light elements and then discuss the 
current status of the observations in relation to theory of each of the light elements. 
Finally, we will discuss the current limits on physics beyond the standard model. 



2 Big Bang Nucleosynthesis 



2. A Brief Primer on the Theoretical Prediction 

Conditions for the synthesis of the light elements were attained in the early Universe 
at temperatures T <, 1 MeV, corresponding to an age of about 1 second. At 
somewhat higher temperatures, weak interaction rates were in equilibrium, thus 
fixing the ratio of number densities of neutrons to protons. At T 3> 1 MeV, (n/p) — 
1. As the temperature fell and approached the point where the weak interaction 
rates were no longer fast enough to maintain equilibrium, the neutron to proton 
ratio was given approximately by the Boltzmann factor, (n/p) ~ e -Am / T , where 
Am is the neutron-proton mass difference. The final abundance of 4 He is very 
sensitive to the (n/p) ratio. 

The nucleosynthesis chain begins with the formation of deuterium through the 
process, p + n — > D + 7. However, because the large number of photons relative 
to nucleons, r\~ x = n^/ns ~ 10 10 , deuterium production is delayed past the point 
where the temperature has fallen below the deuterium binding energy, Eb = 2.2 
MeV (the average photon energy in a blackbody is E 1 ~ 2. IT). When the quantity 
?y _1 exp(— Eb/T) ~ 1 the rate for deuterium destruction (D + 7 — > p + n) finally 
falls below the deuterium production rate and the nuclear chain begins at a tem- 
perature T ~ O.lMeV. 

The dominant product of big bang nucleosynthesis is 4 He resulting in an abun- 
dance of close to 25% by mass. This quantity is easily estimated by counting the 
number of neutrons present when nucleosynthesis begins. When the weak interac- 
tion rates freeze-out, at T « 0.8 MeV, the neutron to proton ratio is about 1/6. 
When free neutron decays are taken into account prior deuterium formation, the 
ratio drops to (n/p) « 1/7. Then simple counting yields a 4 He mass fraction 

[1 + (n/p)] 

In the standard model, the 4 He mass fraction depends primarily on the baryon 
to photon ratio, 7/ as it is this quantity which determines the onset of nucleosynthesis 
via deuterium production. But because the (n/p) ratio is only weakly dependent on 
77, the 4 He mass fraction is relatively flat as a function of 77. The change due to the 
uncertainty in the neutron half-life is small (this effect is shown in Fig. 1). When 
we go beyond the standard model, the He abundance is very sensitive to changes in 
the expansion rate which can be related to the effective number of neutrino flavors 
as will be discussed below. Lesser amounts of the other light elements are produced: 
D and 3 He at the level of about 10~ 5 by number, and 7 Li at the level of 10 -10 by 
number. 

The resulting abundances of the light elements are shown in Figure 1 from the 
calculations in ref. 2. The curves for the 4 He mass fraction, Y, bracket the computed 
range based on the uncertainty of the neutron mean-life which has been taken asa 
r n = 887 ± 2 s. The 4 He curves have been adjusted according to the corrections 
in ref. 4. Uncertainties in the produced 7 Li abundances have been adopted from 
the results in ref. 5. Uncertainties in D and 3 He production are negligible on the 
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Fig. 1. The light element abundances from big bang nucleosynthesis. 
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scale of this figure. The boxes correspond to the observed abundances and will be 
discussed below. It is clear that as the observational boxes line up on top of each 
other, there is an overall agreement between theory and observations in the range 
1710 = 10 10 7? = 2.8 - 4.5. 

3. The Observations 
3.1. 4 He 

Because helium is produced in stars along with heavier elements, it is necessary 
to look for primordial helium in regions where the stellar processing is minimal, i.e., 
in regions where the abundances of elements such as carbon, nitrogen and oxygen 
are very low. The 4 He abundance in very low metallicity regions is best determined 
from observations in extragalactic HII regions of Hell — > Hel recombination lines. 
There are extensive compilations of observed abundances of 4 He, N, and O, in many 
different galaxies Btm In Figure 2, the 4 He vs. O/H data is shown along with its 
associated linear fit. This is an updated version of the plot from ref. 9, including 
data from ref. 8 (details of this fit are given in the last line of Table 1). 



Fig. 2. The observed abundances of 4 He vs. O/H in extragalactic HII regions along with a linear 
fit to the data. 



In Table 1, various fits to the data and subsets of the data are given. Details 
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Table 1. Linear Fits for Y vs. O/H. 



Set 


# Regions 


r 


x'/dof 


Y P 


10 2 x slope 


1 P 


"All" 


49 


0.56 


0.78 


.234 ± 


003 


1.14 ±0.24 


0.239 


1st cut 


41 


0.51 


0.61 


.232 ± 


003 


1.38 ±0.36 


0.238 


-outliers 


34 


0.45 


0.70 


.232 ± 


003 


1.39 ±0.38 


0.238 


2nd cut 


21 


0.41 


0.64 


.229 ± 


005 


2.37 ± 1.13 


0.238 


-outliers 


19 


0.40 


0.70 


.229 ± 


005 


2.42 ± 1.15 


0.238 


C 


22 


0.35 


0.71 


.232 ± 


003 


1.58 ±0.54 


0.238 


IZ + "All" 


55 


0.53 


0.82 


.236 ± 


002 


1.02 ±0.23 


0.240 


IZ + 1st cut 


47 


0.49 


0.69 


.234 ± 


003 


1.24 ±0.33 


0.239 



concerning the subsets of the data shown (as well as a more complete discussion 
on primordial 4 He) can be found in ref. 9. As one can see there is a considerable 
degree of stability in these fits, leading to a 2 a upper limit of 0.238 -0.240 for the 
primordial abundance of 4 Hc. There is in addition an overall systematic uncertainty 
of about 0.005 in Y p . 

The size of the assumed systematic errors has recently been questionedllj'Ell. 
Several sources for the systematic errors have been discussed and it has been ar- 
gued that the systematic error is infact significantly larger than 0.005. Though the 
cumulative effect from several systematic uncertainties is unclear, linearly adding 
the various systematic errors assumes that they are correlated and will surely result 
in an over-estimate of the true error. On the other hand, the effects of detailed 
radiative transfer may be more important than previously thoughttil in determin- 
ing the 4 He abundance, and pending new results concerning these calculations it 
is premature to estimate the size of the correction to Y, if any. For now, we will 
continue to assume that <T S y S tcmatic = 0.005 though we must bear in mind that the 
true error may be somewhat larger. 

As was noted, from Table 1 it appears that the results for Y p , the intercept of 
the linear fit to the Y vs. O/H data, are quite stable. In order to test the stability 
of the fitf we have randomly sampled 49 points from the 49 observed points from 
refs. 6 and 7 and included in fig. 2. 10,000 separate sets of 49 points were obtained. 
We then performed a least squares linear fit to each set of points. The helium 
fraction for the zero metallicity point for each of the fits is illustrated in fig. 3. The 
fits are strongly peaked around the mean value of .232 ± .003 (.005 at 95% CL) in 
remarkable agreement with the simple linear fit in the first line of Table 1. We can 
conclude that there are apparently no observed points that are strongly influencing 
the fit. We have similarly run tests on the first cut data set of 41 points, with the 
resulting mean value for the intercepts of 0.232 ± 0.003 (0.007 at 95% CL). 

Subsequent to the analysis in ref. 9, concerning primordial 4 He, there has been 
some newer dataQ on 4 He in low metallicity extragalactic HII regions. Here we show 
how these results modify the above fits of 4 He vs O/H. The data in Izatov et alfl 
consists of 4 He, O/H, and N/H measurements in 10 extragalactic HII regions. Four 

"Work done in collaboration with G. Steigman 
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Fig. 3. The distribution of intercepts from 10 4 linear fits of Y vs. O/H based on a random selection 
of data from the full data set of 49 extragalactic HII regions used in ref. 9. 

of these regions had been observed in Skillman et al. @ and these have been averaged 
together with the new data. Thus there are a total of 6 new points in the data set. 
We will not enter into the discussion regarding the methods of data reductions used 
by Izatov. Here we have used their data based on having taken electron densities 
from the ratio of collisionally excited SII lines (see ref. 8, for details). Though 
there is a slight shift upwards, the shift is entirely consistent within the statistical 
uncertainties. One could however make a case that based on this set of data (the 
combined set) that the 2 a upper limit to Y p is 0.239 which moves up to 0.244 when 
a systematic uncertainty of 0.005 is added. 

It is also of interest to compare the correlation of 4 He with N/H. Indeed, a com- 
parison of N/H vs O/H is also of interestoij. Since the stars which produce oxygen 
do not coincide exactly with those which produce nitrogen, the linear regression of 
Y vs N/H is another consistency check. For the same data as discussed above, the 
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Fig. 4. The observed abundances of 4 He vs. N/H in extragalactic HII regions along with a linear 
fit to the data. 

corresponding results to Y vs N/H are given in Table 2. Again, we see that there is 
a considerable degree of consistency within the data. The 2 a upper limit is slightly 
higher. The stability of the Y vs. N/H fits was also checked ( by the same method 
used to generate Fig. 3 ). For the full data set of 49 points from ref. 9 (line 1 of 
Table 2), the mean value of intercepts found was 0.234 ±0.002 (0.005 at 95% CL). 
For the first cut set of points (line 2 of Table 2) , the resulting mean value was also 
0.234 with a slightly larger uncertainty ±0.003 (0.006 at 95%CL). 

Looking at Tables 1 and 2, it is not unreasonable to take as the best estimate 
to the primordial mass fraction of 4 He 

Y p = 0.234 ± 0.003 ± 0.005 (2) 

with a 2 a upper limit of 0.240 and 0.245 when systematic uncertainties are included. 
Thus we will assume a range (2 a plus systematic) of 0.223 - 0.245 for Y p which is 
shown in Fig. 1 as the large box bracketing the 4 He curves. 

Before concluding the discussion on 4 He, we would like to upda te an older 
analysiso concerning the correlation of 4 He with C/H. It was suggested^ that be- 
cause carbon (and nitrogen) are produced in intermediate mass stars, that the 
correlation between these elements and 4 He might be more significant than that of 
4 He and O/H. The correlation between 4 He and C/H was tested in Walker et alfl 
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Table 2. Linear Fits for Y vs. N/H. 



Set 


# Regions 


r 


X'/dof 


Y P 


10 1 x slope 


1 P 


All 


49 


0.66 


0.66 


.236 ± 


002 


1.72 ±0.33 


0.240 


1st cut 


41 


0.57 


0.58 


.234 ± 


002 


2.71 ±0.68 


0.239 


-outliers 


34 


0.48 


0.69 


.234 ± 


003 


2.77 ±0.76 


0.239 


2nd cut 


21 


0.47 


0.63 


.231 ± 


004 


4.85 ±2.27 


0.239 


-outliers 


19 


0.44 


0.70 


.232 ± 


004 


4.79 ±2.29 


0.239 


C 


22 


0.46 


0.60 


.233 ± 


003 


3.62 ± 1.17 


0.238 


IZ + "All" 


55 


0.61 


0.74 


.239 ± 


002 


1.51 ± 0.31 


0.242 


IZ + 1st cut 


47 


0.52 


0.70 


.237 ± 


002 


2.35 ±0.64 


0.241 



where the fit to Y vs. C/H gave similar results as for Y vs O/H and N/H albeit with 
larger uncertainties to due to the poor statistics (there were only 6 extragalactic 
HII regions with measured C/H). Recently there have been new measurementstj of 
C/O ratios in seven regions. The data are somewhat better than previous data (at 
least for small C/H) and are shown in Fig. 5. Once again, Y vs C/H is of limited 
value because of the paucity of data, yet yields a consistent (with O/H and N/H) 
fit with an intercept Y p = 0.231 ± 0.004 with a 2 a upper limit of 0.240; the same 
as we have concluded above from O/H and N/H. 



Fig. 5. The observed abundances of 4 He vs. C/H in extragalactic HII regions along with a linear 
fit to the data. 
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3.2. D and 3 Hc 

It is more difficult to compare the primordial deuterium and 3 He abundances 
with the observations. Despite the fact that all observed deuterium is primordial, 
deuterium is destroyed in stars. A comparison between the predictions of the stan- 
dard model and observed solar and interstellar values of deuterium must be made 
in conjunction with models of galactic chemical evolution. The problem concerning 
3 He is even more difficult. Not only is primordial 3 He destroyed in stars but it is 
very likely that low mass stars are net producers of 3 He. Thus the comparison be- 
tween theory and observations is complicated not only by our lack of understanding 
regarding chemical evolution but also by the uncertainties of the production of 3 He 
in stars. Furthermore, due to the large uncertainty in the age of the galaxy, it is 
difficult to directly compare the observed and modeled values of deuterium and 3 He 
since the exact galactic age at which the pre-solar epoch occurs is unknown. We 
will return to this problem below. 

It appears that D/H has decreased over the age of the galaxy. The (pre)-solar 
system abundance of deuterium is inferred from a variety of 3 He measurements in 
the solar wind, lunar soil samples and meteoritic samples. There are two compo- 
nents of 3 He: one associated with the pre-solar abundance of deuterium and 3 He 
(the deuterium was converted in the early premain-sequence stage of the sun); while 
the second is 3 He trapped in meteorites which represents the true solar value of 3 Hc 

= (4.1 ± 1.0) x 1(T 5 (3) 



= (1.5 ±0.3) x 1CT 5 (4) 



The difference_between the two 3 He measurements represents the presolar D abun- 
dance and isN D/H « (2.6 ± 1.0) x 10~ 5 . On the other hand, the present ISM 
abundance of D/H islla 

D/H = 1.60 ± 0.09l;^ x 10 -5 (5) 

It is this lower limit on D/H (since D is only destroyed) that provides us with the 
the upper bound on rj. It is shown as the lower side of the D and 3 He box in 
Fig. 1. Thus, if 7710 is in the range 2.8 - 3.9 then the primordial abundance of 
D/H is between 4.5 - 8 xlO -5 , and it would appear that significant destruction of 
deuterium is necessary. 

Note that there are reported detections of D in high redshift, low metallicity 
quasar absorption systemsEZHia. In principle, these measured abundances should 
represent the primordial one. Unfortunately, at this time, the two observations do 
not yield a single value for D p . In oneEJ, D/H « 1.9 — 2.5 x 10~ 4 while the second 
givesta, D/H « 1 — 2 x 10~ 5 . (One should note that it has been argued that these 
types of measurements are subject to large systematic correctionsEJ, thus it may not 
be all that surprising that the two measurements differ.) The former observation 



D + 3 Hc 
H 
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is shown in Fig. 1 by the small box on the D/H curve at a value of 7710 « 1.5. As 
one can see the corresponding value of Y p (at the same value of 77) is in excellent 
agreement with the data. 7 Li is also acceptable at this value as well. Due to the 
still some what preliminary status of this observation it is premature to fix the 
primordial abundance with this value. A high value for the D abundance would 
require an even greater degree of D destruction over the age of the galaxy. The 
lower measurement for D/H is problematic for both 4 He and 7 Li. 

To trace the evolution of D and 3 He, we can make use of some simple models 
of galactic chemical evolution. For a general review on chemical evolution see ref. 
20. We shall adapt her same basic notation in the following discussion. One can 
relate the deuterium abundance at any time to the gas mass fraction noting that 
deuterium is only destroyed in stars. The equation for the gas mass in the galactic 
disk is 

^E=e(t)-m + f(t)- (6) 

In this equation, e(t) is the gas ejected from stars at the end of their lifetimes either 
in planetary nebulae or supernova events, 

TTL U pp 

e(t) = J (m — to.r)V>(* — T m )<f>(m) dm. (7) 

m(t) 

In these equations, f(t) is the infall rate onto the disk of gas whose origins is external 
to the galaxy. </>(m) is the initial mass function (IMF), and ip(t) is the star formation 
rate (SFR). m up p is the upper mass limit on <f){m). r m is the lifetime of a star of 
mass to, and is the remnant mass left over by a star of mass to. m(t) is is the 
mass of a star which at time t is returning gas back into the ISM. 

In order to model the deuterium evolution, we consider D to be totally astrated 
within stars. We may extend equation @, 

^=_e(i)D + /(t)(D / -D). (8) 

where Df is the mass fraction of D in the infall gas. Using the instantaneous 
recycling approximation (an assumption which ignores delays due to finite stellar 
lifetimes, and is an acceptable approximation at early times and for elements which 
are produced in massive stars whose lifetimes are short) we can combine Eqs. (Q) 
and (|^) to give the total amount of deuterium destruction relative to its initial 
primordial value dJ3 

D = a R/(l-R) (9) 
Dp 

where a is the gas mass fraction and the return fraction R is given by 



R= J (m — mn)4>(m) dm (10) 

m(t) 
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By a suitable choice for lhe_IMF (and SFR) a considerable amount of deuterium 
may have been destroyedE3'Ej. 

As we noted above, to match the abundances in a model of galactic evolution 
to the observations requires knowledge of the time at which the solar system was 
formed. In order to solve this problem, wet-3 have incorporated the production of 
the long-lived r-process nuclear chronometers 232 Th, 235 U, and 238 [/ into chem- 
ical evolution models. Previously, these chronometers have been used in simple 
chemical evolution models in an attempt to constrain the over-all galactic ageEl 
This is accomplished by finding a time at which the modeled ratios 232 Th/ 238 U 
and 235 ^/ / / 238 [/ match (if ever) the values of these ratios determined from meteorite 
data { 232 Th/ 23S U = 2.32 and ^U/ 23S U = .3170). This time is taken to be the 
pre-solar epoch, and 4.6 Gyrs. later, today. 

A full description of how nuclear chronometers may be incorporated in general 
into galactic chemical evolution models can be found in ref. 25. We summarize the 
main points below. Nuclear chronometers are incorporated into the calculation by 
an extension of the basic equation of the gas mass in the galactic disk discussed 
above. We may extend equation (^) to determine the rate of change in the number 
of nuclear species A, 

^ = ™-^ + ^ + 5 £>-^ ■ (U, 

In this equation, Pa is the number of newly synthesized nuclei of species A per 
unit mass going into star formation. The relative production ratios, -P232/-P238J and 
-F235/-P238 > have been determined by supernova modelsc3, P232/-P238 = 1-60, and 
-P235/-P238 = 1-16. Since we have assumed the infall gas to be primordial, Zf, the 
metallicity of the infall gas, has been taken to be zero. And finally, Xa is the rate 
of decay of nuclear species A. As a simplification, the amount of nuclear species 
"A" which decays while locked up in stars was not taken into account equation (3). 
Only large stars will contribute to the production of these nuclear species and they 
have lifetimes which are very short when compared with the decay times. 

Equations (0) and ( |Ti"| ) can then be solved numerically. In order to satisfy obser- 
vational constraints, the value of Na(0) is allowed to be non-zero. This corresponds 
to an initial enrichment (So) in the galactic disk of metals. This solves the so-called 
G-dwarf problem; where observationally, there appear to be very few metal poor 
dwarf stars in the diskcH. 

We have found a wide range of chemical evolution models capable of destroying 
D by a factor of 5, which also satisfy our nuclear chronometer constraint. These 
models result in a range of galactic ages (11.1 — 15.0 Gyrs). The result of one of 
these models is shown in Fig. 6. 

There are however potential problems for 3 He. A lower limit on rj was derivedEa 
by noting that although stars can destroy 3 He, even very massive stars, still preserve 
at least 25 % of the initial D + 3 He. (It is the sum of D and 3 He that is important as 
D is burned to 3 He in the premain-sequence phase of stars.) A value of 7710 lower than 
2.8 would yield (D+ 3 He) /H > 10 -4 so that even if the maximal amount of 3 He is 
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Fig. 6. The evolution of D/H in several models of galactic chemical evolutio: 




destroyed, it would still exceed the presolar valued of (D+ 3 He)/H w 4.1±1.0x lflr 5 . 
But, in low mass stars, 3 He is produced rather than destroyed. A calculatiorO of 
the final 3 He abundance relative to the initial abundance of D + 3 He gives, 

( 3 He/H)/ = 1.8 x 10- 4 (^f ) + 0.7 [(D + 3 He)/H] . (12) 

so that a 1 M Q star produces 2.7 times as much 3 He as the initial D -t-. 3 He (for (D 
+ 3 He)initiai = 9 x 10~ 5 ). This would lead to an evolutionary behavioicj of the type 
shown in Fig. 7. The chemical evolution model has been chosen so that D/H agrees 
with the data and assumes that r) 10 = 3. The problem being emphasized concerns 
3 He and can be seen by comparing the solid curve with the filled diamonds. 

A number of questions regarding the 3 He discrepancy can be raised. First, 
one can ask whether or not the formula for 3 He production is valid. There is 
indeed evidence that 3 He is produced in planetary nebulaeu where abundances 
are found as high as 3 He/H ~ 10 . Secondly, one may ask how uniform are the 
3 He measurements. As it turns out, they are in fact not very uniformES and show 
variations as large as a factor of 5 between different HII regions. Indeed, there may 
even be acorrelation between the size (mass) of the region and the amount of 3 He 
observedti3. Thus it may be possible that the solar values are depleted in 3 He. It is 
hoped that future 3 He observations will help to resolve this puzzle. 
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Fig. 7. The evolution of D/H (dashed curve), 3 He/H (solid curve) and (D + 3 He)/H (dotted 
curve) as a function of time. Also shown are the data at the solar epoch t Rj 9.6r£lyr and today 
for D/H (open squares), 3 Hc/H (filled diamonds) and (D + 3 He)/H (open circle) ,PH . 



Before we conclude with D and 3 He, we address the overall issue regarding 
consistency of the data and the big bang nucleosynthesis calculations. Recently, it 
has been suggestedE 3 ] that there is a crises in big bang nucleosynthesis because of 
an underlying inconsistency between the "predicted" value of rj (from deuterium) 
and the 4 He observations. The conclusion was that to a high degree of confidence, 
the standard model with three light neutrino flavors is excluded. Though we will 
return to neutrino limits in the next section, we can here dispell the notion of a major 
inconsistency. The crises in ref. 33, arises from using a pedagogical expressiorO@ 
relating deuterium and 3 He at a given time assuming that a fraction / of the total 
gas mass passed through stars at time t, and that a fraction g 3 of the initial D + 
3 He survives as 3 He when 3 He production is neglected 

<m + i(!S£u fi^vm da) 



h w, » v h a v 93 r V H , , 

This expression has been usedEl extensively to place a lower bound on rj by inserting 
the observed values for D/H and 3 He/H at the solar epoch. This is what leads to 
7710 > 2.8 when g 3 > 0.25 and corresponds to the upper limit on the D + 3 He box 
in Fig. 1. To be sure, this expression indicates that the best fit! to deuterium is a 
low value for D/H (w 2.3 x 10~ 5 ). Clearly for the corresponding value of 77, there is 
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a significant discrepancy with 4 He. However, rather than to conclude on this basis 
that there is a problem with the standard model, or that we have learned something 
about neutrino physics, it is evident to us that that what we have learned is that 
there is something wrong with the expression (|l3|). Indeed it is known that even 
simple models of galactic evolution with 33 ~ 0.3 (Le. neglecting 3 He production) 
can give a good match to the D and 3 He solar datao. In this case, the discrepancy 
between the prediction and observations of 4 He is at most at the 2 a level. Of course 
when 3 He production is included, the problem is more acuteEj requiri ng , perhaps, 
more dramatic changes in the chemical evolution models of the galaxyO. 

3.3. 7 Li 

Finally we turn to 7 Li. Over the last several Years, there has been a consid- 
erable increase in the number of 7 Li observations c3. 7 Li in old, hot, population 
II stars, is found to have a very nearly uniform abundance. For stars with a sur- 
face temperature T > 5500 K and a metallicity less than about l/20th solar, the 
7 Li abundance shows little or no dispersion beyond what is consistent with errors 
of individual me asure ments. A large quantity of data has been obtained from a 
variety of sourcescSEU and consists of observations of 55 different halo stars. The 
corresponding mean 7 Li abundance is [Li] = 12 + log Li/H = 2.08 ± 0.02 or the 
abundance by number is Li/H — (1.22 ± .04) x 10~ 10 . Systematic errors however 
dominate the uncertainty in the 7 Li abundance. 

Furthermore there has been considerable attention to trends in the data. Namely, 
it appears that there may bejieal correlations in the 7 Li abundances with respect to 
temperature and metallicityEUH. It has been argued that to avoid any temperature 
dependence which may mask a real metallicity dependence, one should renormalize 
the derived 7 Li abundances to a single temperature, say 6200K. For the data dis- 
cussed above, there is only a weak dependence on temperature, [Li] = 1.20 ± 0.44 
+ (14 ±7) x 10~ 5 T. This fit is a modest improvement (at the 87 % CL) over the 
straight average. "Correcting" the data by sliding each point up to 6200K, along 
the line given by the fit above, one finds on the average a slightly higher 7 Li abun- 
dance, [Li] = 2.13 ± 0.02. If one were to take the raw data, (ie. 7 Li abundances 
uncorrected for differing temperatures) then only a questionable correlation with 
metallicity would result [Li] = 2.14 ± 0.05 + (0.025 ± 0.023)[Fe/H], where [Fe/H] is 
the log of the iron abundance relative to the solar iron abundance. In contrast, when 
the corrected lithium abundances are used a more significant correlation appears, 
[Li] = 2.24 ± 0.05 + (0.051 ± 0.022)[Fe/H]. This fit represents an improvement at 
the 98% CL over the straight average. 

The derived 7 Li abundance is quite sensitive to the assumed temperatures of the 
stars. Most authors have obtained abundances through color-based temperatures 
from photometry. In addition, the results depend on the particular stellar model 
used. In Thorburn'sEa analysis of over 70 halo stars, she finds a systematically 
higher temperature and hence 7 Li abundance. In this sample [Li] = 2.26 ±0.01 or 
Li/H = (1.82 ± .04) x 10~ 10 . 
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Fig. 8. 7 Li data shown as a function of temperature. [Li] = 12 + log Li/ H. The horizontal line is 
the mean of this data. 

The same correlation (and to higher degree) between Li and temperature and 
metallicity is present in the Thorburn data. The correlation with temperature in 
that data appears to be very strong [Li] = 0.06 ± 0.37 + (36 ±6) x 10~ 5 T which 
leads to a higher mean 7 Li abundance [Li] = 2.33 ± 0.01. There is also marked 
effect on the fit with respect to [Fe/H]. Using the raw 7 Li abundances, [Li] = 2.49 
± 0.09 + (0.087 ± 0.036) [Fe/H] and using the corrected 7 Li abundances gives [Li] 
= 2.73 ± 0.08 + (0.16 ± 0.03)[Fe/H]. These fits all have significantly lower x 2 's 
than does the straight average. Of course these trends make it more difficult to 
determine with any certainty the true primordial abundance. 

Recently however, the method oLobtaining temperatures by means of broad- 
band photometry has been criticized^ and a strong case has been made to usa_a 
spectroscopic method based on Balmer line profiles. This method was employedllj 
in reanalyzing a relatively large set of halo stars (24 with T > 5700K with [Fe/H] 
< -1.4). In this data set the mean value for 7 Li is [Li] = 2.21 ± 0.02 or Li/H 
= (1.62 ±0.07) x 10~ 10 . However in this case, there is absolutely no evidence 
for a trend with respect to temperature, [Li] = 1.90 ± 0.57 + (5 ±9) x 10~ 5 T 
which corrects the mean value to [Li] = 2.22 ± 0.02. More importantly there is no 
correlation with metallicity either. The uncorrected data gives [Li] = 2.20 ± 0.10 
+ (0.00 ± 0.04) [Fe/H] while the "corrected" data gives [Li] = 2.23 ± 0.10 + (0.01 
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± 0.04) [Fe/H]. 

Though it is difficult to place a specific value on the total systematic error in 
any of the 7 Li abundances, from the spread in the obtained values it would appear 
that a fair estimate of the systematic uncertainty is about 0.1 dex. It may indeed 
be smaller in the last analysesEl. Taken together with what should be considered 
the best estimate for the 7 Li abundance in halo dwarfs, we arrive at a final 7 Li 
abundance of 

[Li] = 2.21 ±0.02 ±0.1 (14) 

corresponding to a range for Li/H between ~ 1.2 and 2.2 x 10~ 10 which is shown 
in Fig. 1 as the observed range for 7 Li. These values (and their uniformity) should 
be compared with observations of Li in younger stars where the abundance can be 
much larger (by an order of magnitude) and shows considerable dispersion. Two 
key questions remain however: how much of the observed Li is primordial (since Li 
is known to be produced) ; and how much of the primordial Li remains in the stars 
where Li is observed? 

Aside from the big bang, Li is produced together with Be and B in cosmic ray 
spallation of C,N,0 by protons and a-particles. Li is also produced by a — a fusion. 
Be and B have recently been observed in these same pop II stars and in particular 
there are a dozen or so, stars in which both Be and Li have been observed. Thus 
Be (and B though there is still a paucity of data) can be used as a consistency 
check on primordial Li c3. Based on the Be abundance found in these stars, one can 
conclude that no more than 10-20 % of the 7 Li is due to cosmic ray nucleosynthesis 
leaving the remainder (an abundance near 10 -10 ) as primordial. It is also possible 
however, that some of the initial Li in these stars has been depleted. Standard 
stellar modelsEJ predict that any depletion of 7 Li would be accompanied by a very 
severe depletion of 6 Li. Until recently, 6 Li had never been observed in hot pop II 
stars. The observationEl of 6 Li (which turns out to be consistent with its origin in 
cosmic-ray nucleosynthesis and with a small amount of depletion as expected from 
standard stellar models) is a good indication that 7 Li has not been destroyed in 
these starsc3. 

3.4. Summary 

Consistency of the standard model of big bang nucleosynthesis relies on the 
concordance between theory and observation of the light element abundances for a 
single value of rj. We now summarize the constraints on 77 from each of the light 
elements. From the 4 He mass fraction, Y < 0.240(0.245), we have that 7710 < 
2.9(4.5) as a 2a upper limit (the higher value takes into account possible systematic 
errors). It is appropriate to note here the sensitivity of Lhe bound on 77 to the 
assumed 4 He upper limit. When the data of Izatov et alJj are not included, our 
upper limit to Y waso Y < 0.238(0.243). In this case the corresponding upper 
limit to 77 is 7710 < 2.5(3.9). Because of the sensitivity to the assumed upper limit 
on Y p , the upper limit on 77 from D/H, though weaker is still of value. From D/H 
> 1.5 x 10 -5 we have 7710 < 7. The lower limit on 77, comes from the upper limit on D 
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+ 3 He and is 7710 £ 2.8 if one ignores 3 Hc production. We stress that this limit is only 
exemplary as the upper limit on D + 3 Hc depends critically on models of galactic 
chemical evolution, which are far from being understood. Finally, 7 Li allows a broad 
range for 77 consistent with other light elements. When both uncertainties in the 
reaction rates and systematic uncertainties in the observed abundances are taken 
into account, 7 Li allows values of 7710 between 1.3 - 4.9. The resulting consistent 
range for 7710 becomes 2.8 - 4.5. These bounds on 77 constrain the fraction of critical 
density in baryons,f2s to be 

0.01<n B <.l (15) 

for a hubble parameter, h Q , between 0.4 - 1.0. The corresponding range for ilsh 2 , 
is 0.010 - 0.016. 

4. Constraints on Physic beyond the Standard Model 

Limits on particle physics beyond the standard model are mostly sensitive to the 
bounds imposed on the 4 He abundance. As is well known, the 4 He abundance is pre- 
dominantly determined by the neutron-to-proton ratio just prior to nucleosynthesis 
and is easily estimated assuming that all neutrons are incorporated into 4 He (see 
Eq. (P). As discussed earlier, the neutron-to-proton ratio is fixed by its equilibrium 
value at the freeze-out of the weak interaction rates at a temperature Tf ~ 1 MeV 
modulo the occasional free neutron decay. Furthermore, freeze-out is determined 
by the competition between the weak interaction rates and the expansion rate of 
the Universe 

G F 2 T f 5 ~ r wk (T» = H(Tf) ~ y/G N NT f 2 (16) 

where N counts the total (equivalent) number of relativistic particle species. The 
presence of additional neutrino flavors (or any other relativistic species) at the time 
of nucleosynthesis increases the overall energy density of the Universe and hence the 
expansion rate leading to a larger value of Tf, (n/p), and ultimately Y p . Because of 
the form of Eq. (|l6|) it is clear that just as one can place limits c3 on N, any changes 
in the weak or gravitational coupling constants can be similarly constrained (for a 
recent discussion see ref. 46). 

In the standard model, the number of particle species entering into eq. ( |l6| ) can 
be written as N = 5.5 + %N V (5.5 accounts for photons and e^). The observationally 
derived primordial 4 He abundance from Eq. of Y p = 0.234 ± 0.003 ± .005 
translates into a best value for N v = 2.2 ± 0.27 ± .42 (assuming a central value for 
7710 of 3.6, with 7710 = 3, the best value for N v increase to 2.4) which implies a 2a 
upper limit of 2.74 which is extended to N v < 3.16 when systematics are included. 
At face value, such a limit would exclude even a single additional scalar degree pi 
freedom (which counts as |) such as a majoron unless it decoupled early enoughcJ 
so that its temperature, Tg at the time of nucleosynthesis was suppressed so that 
{Tb/T v ) a < j(-16) = -28. In models with right-handed interactions, and three right- 
handed neutrinos, the constraint is more severe. The right-handed states must have 
decoupled early enough to ensure {T UR /T UL ) A: < (.16)/3 ~ .05. The temperature 
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of a decoupled state is easily determined from entropy conservation, (T x /T v ) — 
((43/4)/iV(T ti )) 1/3 . One additional scalar requires N(Td) > 28 or decoupling must 
have occurred above the QCD phase transition at a temperature Td > Xqcd ~ 200 
MeV. Three right-handed neutrinos would require N(Td) >, 100, so that Td > 
Mw If right-handed interactions are mediated by additional gauge interactions, 
the associated mass scale becomes Mz> > 0(100) TeV! 

The limits on N v , however, are sensitive to the upper limit on 4 He which is 
in turn sensitive to assumed systematic errors and to the lower bound on rj. In 
addition, the limits described above may be overly restrictive^. The best value for 
N v is 2.2 and may in fact be unphysical if v T is lighter than ~ 1 MeV, as is quite 
likely. In this case the limits on N u must be accordingly renormalized as has been 
discussed extensively in ref. 48. In Fig. 9, the effect of renormalizing the limit on 
N v is shown. In Fig. 9, the central value of Y p — 0.232 was assumed. For the higher 
value of 0.234, an additional « 0.2 should be added to the limits on N v . 



Fig. 9. The 95 % CL upper limit on N v as a function of the systematic uncertainty in Y p . The 
solid (dashed) curves correspond to the condition that 2V„ > 2 (3). Each of these two cases is 
shown for three choices of r?io: 1.5, 2.8, and 4.0. The smaller values p£?7 yield weaker upper limits 
and systematic errors have been described by a top-hat distributiorcj ■ 



5. Summary 

In summary, we have argued for the overall agreement between theory and ob- 
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servations as they pertain to the light element abundances as well as concordance 
between big bang nucleosynthesis and galactic cosmic-ray nucleosynthesis. There is 
now quite a bit of data on 4 He which appears to be coherent. Because the inferred 
primordial 4 He abundance found from extrapolating the data to zero metallicity 
gives a somewhat low value, smaller values for the baryon to photon ratio are pre- 
ferred. In contrast, there may be a problem with 3 He. The data in the ISM shows a 
large dispersion, which may be an indication of possible local pollution, and spans 
a factor of ~ 5 in abundance. Additional data is clearly needed here before strong 
conclusions can be made regarding potential problems. The solar value of 3 He 
seems to indicate that little 3 He is produced over the early history of the galaxy. 
While higher values of the baryon to photon ratio certainly help in this direction, 
we believe that this is ultimately a problem to be resolved by chemical evolution 
models. 7 Li, because of relatively large allowed range (dominated by systematic 
uncertainties) allows a wide range for the baryon to photon ratio. 

There are however, open issues: Are the quasar line-of-sight measurements0lll 
of D/H real; Why isn't there more 3 He, particularly in the solar system; Can the 
statistical and systematic errors in He measurements be reduced; Can the large sys- 
tematic errors in the 7 Li abundance be reduced? Clearly new data will be necessary 
to resolve these problems. Nevertheless, in spite of these uncertainties, nucleosyn- 
thesis is still able to set strong constraints on physics beyond the standard model. 
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